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CROP STATUS EVALUATIONS AND YIELD PREDICTIONS 


Principal Investigator: J. R, Haun 
Department of Horticulture 
Clemson University, Clemson, S. C. 

Introduction - The work reported herein was based on the program outlined in the reference 
NASA Headquarter's Y-No, 41-001-043. During the course of this work it was decided 
to extend the contract to include the remaining portion of the 1974-75 winter wheat 
season and to collect data on a second season of both spring and winter wheat. 

Administrative and personnel arrangements - Interviews prior to employment and subsequent 
training procedures were carried out for the following positions reporting to the Principal 
Investigator at Clemson. 

Research Analyst II - The person in this position has a BA degree with major 
in mathematics. 

Graduate Assistants - One, with a BS degree is a PhD candidate and has been 
active during the report period. Another, with an MS degree, also a PhD 
candidate, has been tentatively accepted to work on this contract. 

Laboratory Technician - University students have been employed part time 
during the school term and full time during the summer months. 

Data collection - Arrangements were made with the following stations of the University 
of North Dakota: Minot, Dickinson and Williston, and with an ASCS employee at 
Williston to collect daily data on spring whea* development and environmental variables. 
Spring wheat plots were also grown In the field and in the greenhouse at Clemson to 
provide data from conditions widely different climatically from commercial spring wheat 
areas, hor winter wheat data, arrangements were made with the following stations of 
Kansas State University: Colby, Garden City, Parsons and Manhattan. Winter wheat 
plots were also planted in November 1974 at Clemson and the collection of growth data 
is currently still being taken for the 1974-75 season at all locations. Visits were made 
to all locations 'O discuss procedures, observe experimental conditions, and provide 
instructions for ofcse.vers. 

Tables 1-5 of the Appendix present the data collected on spring wheat for the first 
season (1974). The variable GR0 represents the daily growth rate as calculated for the 
average change in growth stage from the previous day of 25 tagged plants. The remaining 
abbreviations ore explained in the list preceding these tables. 
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Long-term weather records for North Dakota and Kansas were obtained on magnetic tapu 
from the National Climatic Center at Asheville, N. Carolina. Much work was necessary 
to develop procedures for decoding and/or interpreting these records. 
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Plan of data processing and analytes - Figs. 1 and 2 detail the various steps or phases of 
programming and analysis, with an indication of their sequence. At the date of this 
report items 1-15, and 19 for 1974 in North Dakota (Fig. 1 is for spring wheat analyses 
which incljde the data from Clemson) and items 1, 2, 4, 5, 7, and 8 for 1974-75 Kansas 
(Fig. 2 is for winter wheat also including Clemson data) have been obtained and/or 
completed. 

Progress and results on specific phases - 

A/ Comparison of growth-environment relationships for wheat, barley and rye - Figs. 3 
and 4 show the comparative growth rates of these crops in greenhouse and field conditions. 
Obviously differences in response among these cereals are very small. 

Growth rate prediction model- were developed for the greenhouse data on wheat, barley 
and rye by the step-down multiple regression procedure. This procedure begins with all 
of the independent variables and sequentially removes variables that are below a prede- 
termined level of significance after each step. The daily growth rate (see Fig. 3) was 
used as the dependent variable and the following basic independent variables were used: 
maximum temperature (TX), minimum temperature (TN), Langleys (LN), relative 
humidity hours above 30% (RH), soil temperature at 8:00 a.m. (ST), degree hours above 
15.6° C. (60° F. ) (DH). The data for each basic variable for each of the two preceding 
days were also used as independent variables (i.e. lagged variables). The basic and 
lagged variables were used as additional independent variables in the form of the following 
transformations: square and cube; and cross products of all basic and lagged variables 
that were relatively more prominent in preliminary analyses using only untransiormed, 
squared, and cubed variables. Cumulative growth (CUM) was used as an independent 
variable in the total of 133 variables. 

Tables 1,2, and 3 show the statistics for the final steps of these analyses. These three 
analyses are similar in that they contain the same number of independent variables, 
similar F and values, and are based on the same number of degrees of freedom. They 
are dissimilar with respect to the nature of the independent variable'.. However, since 
many of these variables closely resemble each other (e.g. LN12 anti LN13) this is not 
conclusive evidence that these combinations represent significant differences in the 
response of these crops to the ’iwironment. 

Growth rate prediction models were developed for Clemson field data of 1971-72 for 
wheat, barley and rye by the step-down multiple regression procedure. The daily growth 
rate was used as the dependent variable (see Fig. 4) and the same independent variables 
as for the greenhouse study except soil temperature. Lagged and transformed variables 
were also used as in the greenhouse study. 

Tables 4, 5, and 6 show the statistics for the final steps of these analyses. An evaluation 
of these results with tests of prediction equations on the other two crops are presented in 
Figs. 8, 9, and 10. Figures 5, 6, and 7 show the observed and predicted growth rates 
for the equations in Tables 4, 5, and 6. These data are also presented in the form of dot 
diagrams in the upper left portion of Fig. 8, upper right portion of Fig. 9, and lower 
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left portion of Fig. 10 with the correlation coefficients shown for the relationship of 
predicted to observed, and the regression line for predicted on observed. The remaining 
portions of Figs. 8, 9, and 10 show results of the equation applied to the other two 
crops. In each Fig. it may be noted that the coefficient is lower when the equation 
for one crop is applied to another, e.g. in Fig. 8, r = .926 for wheat on wheat, .863 
for wheat on barley and .650 for wheat on rye. It may be concluded from this study that 
the responses of these three crops to environmental factors are quite similar. Although 
the equations involve different variables they are not necessarily unique for each crop. 
Additional testing would be necessary to confirm significance of the particular combinations. 
In the event it becomes possible to differentiate acreage of these three crops by remote 
sensing it will be desirable to perfect appropriate equations for each crop. 

B, Ev aluation of plant observation techniques - This study was made on spring wheat 
grown at Clemson to further evaluate^ the established procedure for daily plant observations. 
The usual procedure has been to observe the growth stage (number of leaves and fraction 

of the emerging leaf) on a single culm (stalk), the first to develop, on each of 10 or more 
individual plants. The average of the 10 culms is then used as the stage of development. 
Since each plant later produces many tillers (branches) that also become culms quite 
similar to the original one, but with fewer nodes, the adequacy of a single culm obser- 
vation per plant was questioned, i.e. does the growth stage of a single culm sufficiently 
represent the growth stage of an entire plant with 5 to 10 tillers? In order to answer 
this question, all of the tillers on 10 plants were observed daily for 25 days. Growth 
rates calculated from the 10 main culm observations (one culm per plant) were compared 
with growth rates involving the tillers combined with main culms. A graph of these 
calculations is presented as Fig. 11. In general, the two methods compare favorably. 

The number of tillers is indicated in the lower line of this graph. This number varied 
since some of the tillers died and later others developed. Since the growth rates for the 
two methods are not completely identical another study is currently under way on winter 
wheat. It is possible that the sample size for the main culm method was not large enough 
in the study on spring wheat. All of the growth observations from field locations in 
North Dakota and Kansas (1974) are now based on the main culms of 25 plants per 
location. Based on previous experience this should provide excellent indices of plant 
development. 

C. De velopment of model for predicting the date on which 50% of the spring wheat crop 
is planted - Eight years of ASCS-USDA weekly crop status reports were checked for 
data on percentage of acreage planted in each division. These data were regressed on 
days of the year in order to place regression lines and then select an estimated time of 
50% planting for each division. An example of this procedure is included as Fig. 12. 

The regression line was based on the five reports indicated by the plotted points. The 
estimated time of 50% planted is shown by the circled point on the line corresponding to 
April 27 (day 1 17 of the year). These estimated dates (for 9 divisions for 7 years, leaving 
one year for a test) were then regressed on 49 variables by the step-down procedure. Data 
for 1973 was omitted for use in testing the later model. 1974 data were used in the analysis 
since the year was very atypical and therefore broadened the scope of data for the analysis. 
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The following basic independent variables were used: running three, six, -and nine day 
sums of average minimum and maximum temperature values (C°) (N3S, N6S, N9S, X3S, 
X6S, X9S), estimated soil moisture (EO in %), preseason precipitation (PP in cm), and 
location (LOC). The basic variables were used as additional independent variables in 
the form of the following transformations: square and cube; and cross products of all 
basic variables. Table 7 shows the statistics for the final step of the analysis, with 
Fig. 13 showing the relationship of the predicted values to the actual values. 

The prediction equation (Predicted planting date = 136.7 + 0.0055(E0 X N9S) - etc. ) 
was tested on daily data not used in the analysis for 1973. An example of the predicted 
values for one division (represented by dots) is presented ?n Fig. 14. The curved line 
is the line of best fit (considering only 1st and 2nd order equations) through the pointu. 
The straight line (45°) passes through points where actual and predicted values are equal. 
The objective of the procedure is to find the date where the value predicted by the 
equation is equal to the actual Julian date. Thus, the point where the regression line 
crosses the 45° line is chosen as the predicted planting date. For the Northwest division 
of North Dakota, the actual date was missed by 8 days using this method. Results for all 
divisions are as follows: 


Predicted 


Actual 


Division 

Julian Date 

Julian C 

Northwest 

123 

131 

North Central 

138 

133 

Northeast 

124 

127 

West Central 

116 

124 

Central 

138 

118 

East Central 

124 

111 

Southwest 

126 

126 

South Central 

100 

116 

Southeast 

126 

112 


Avg. 124 

122 


The accuracy of this model is considered satisfactory for the purpose of finding the 
historic dates (i.e. years when no USDA reports were available) on which to initiate 
our yield prediction model. The model was used to obtain divisional dates of planting for 
the years 1965 and 1966. The predicted d visional average Julian day of 50% planted 
was 143 (May 23) for 1965 and 119 (April 2 r , 1966. These values compare with 
133 for 1965 and 125 for 1966 as calculated Horn the recent USDA report, N. D. Wheat 
Historic Estimates 1955-1970, Ag. Statistics No. 33. The values for 1965 and 1966 were 
added to the values for 1^67-72 to obtain 72 divisional planting dates on which to begin 
calculations for the weekly growth indices represented by box 17 in Fig. 1. 1973 and 1974 
data were omitted since the tapes obtained from Asheville do not contain daily weather 
data for these years. 
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D. Development of growth-rate prediction models for spring wheat - Models have been 
dev eToped for spring wheat for weekly periods of the 1974 season "(kox 16 of Fig. 1 ). 

These were developed from analysis of 1974 data from Dickinson, Willlston (2 locations), 
Minot, North Dakota and Clemson, South Carolina. Since observations were not begun 
on the exact date of emergence at each location, It was necessary to establish a stage of 
development In common for all locations on which to base the subsequent weekly Increments 
of data for analysis. The stage 2.4 was used as the point In common to represent the end 

of the first week of development. Successive 7-day increments were used for the subsequent 
weekly an lyses. Each analysis combines the current week's data with all previous weeks. 
Growth rates vere regressed on the following independent variables: maximum air 
Icn.pcioture (C°), minimum air temperature (C°), precipitation (cm), estimated soil 
moisture (100% by Thornth volte method), solar radiation (Langleys); at 0, 1, 2, and 3 
day log periods; squared, cubed; and 45 selected cross products of the basic lagged 
variables; for a total of 93, by the step-up multiple regression procedure. This procedure 
selects the most highly correlated individual independent variable in the first step, and 
in successive steps variables are added or removed to maintain significance at a designated 
level (considering all possible combinations). Due to the fact that errors have been 
detected in the processing of this data, as well as in the initial stages of the analyses, 
all models previously reported must be corrected. Analyses are currently being repeated 
with the necessary changes. 

E. Developm ent of model for predicting the dote when spring wheat emerges in North 
Dakota - Value; representing initial emergence dates were obtained by reading the SRS- 
USDA Crop Status Reports, 1967-74, for North Dakota. When possible, values were 
obtained for each county within a division (some counties never made mention of emergence). 
These values were then averaged within each division to obtain the day of initial emergence. 
The exactness of these values, however, cannot be determined since the terminology used 

by the reporting agents varied greatly and interpretation of these reports is quite subjective. 
Statistics for this model are presented in Table 8, with Fig. 15 showing the relationship of 
the predicted values to the actual values. The model was based on data for 1967-72, 74, 
leaving 1973 as the test year. An example of the predicted values for one division 
(represented by dots) is presented In Fig . 16. The curved line is the line of best fit 
(considering only 1st and 2nd order equations) through the points. The straight line (45°) 
passes through points where actual and predicted values are equal. The objective of the 
procedure is to find the date where the value predicted by the equation is equal to the 
actual Julian date as in the case of the 50% planting date model. Thus, the point where 
the regression line crosses the 45° line is chosen as the predicted emergence dare. The 
test yielded the following divisional results: 



Predicted 

Actual 

Division 

Julian Date 

Julian Date 

Northwest 

153 

137 

North Central 

160 

140 

Northeast 

146 

130 
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West Central 

142 

130 

Central 

152 

119 

East Central 

144 

133 

Southwest 

140 

128 

South Central 

140 

126 

So^.'heast 

141 

130 


Avg. *745" 

W 


Since all of the predictions were too high (average 16 days) It Is probable that some 
standard adjustment will have to be applied if sufficient basis can be obtained. The 
accuracy of this model will be further considered. 
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Table 1. Analysis of variance, regression coefficients, and statistics of fit for the 
dependent variable growth rate of Blueboy II wheat In Clemson greenhouse, January- 
Aprll 1974. 


Source 

df 

MS 

f 

Prob.> F 

— 

Regression 

7 



■|K|1M| 


Error 

39 

0.000631 



■* 

Corrected total 

46 







Partial Regression 
Coefficients 

Student's t 
for HO: B=0 

Probability 

> It| 

Intercept 

0. 15037445 

7.01 

0.0001 

TX2LN1 

0.00001765 

4.61 

0.0001 

LN12 

-0.00000193 

4.95 

0.0001 

RH32 / 1000 

-0.00001435 

2.42 

0.0203 

LN02 / 1000 

0.00040631 

2.40 

0.0215 

TX2ST2 

-0.00015439 

2.50 

0.0169 

DH2 

0.00057213 

1.94 

0.0601 

DH2TN1 

-0.00002263 

1.67 

0. 1034 


The first number following each basic variable indicates the lag period. The second 
number indicates that the variable is squared (?) or cubed (3). Two variables combined 
indicates a cross product. 



Table 2. Analysis of variance, regression coefficients, and statistics of fit for the 
dependent variable growth rate of Keowee barley in Clemson greenhouse, January- 
April 1974, 


Source 


df MS 

F 

Prob.> F 

R2 

Regression 

Error 

Corrected total 

- 

7 0.004246 
39 0.000599 
46 

7.094 

0.0001 

0.560 



Partial Regression 
Coefficients 

Student's t 
for HO: B=0 

Probability 

>|t| 

Intercept 

0.01598093 

0.48 

0.6326 

RH1 

0.00005502 

2.79 

0.0080 

ST 13/ 1000 

0.03579227 

5.64 

0.0001 

ST1TN2 

-0.00061063 

4.21 

0.0001 

LN13/ 1000 

-0.00000378 

4.18 

0.0002 

LN1 

0.00234788 

3.91 

0.0004 

CUMBA 

0.00490153 

3. 15 

0.0031 

ST1LN1 

-0.00012143 

3.12 

0.0034 


The first number following each basic variable indicates the lag period. The second 
number indicates that the variable *-^uared (2) or cubed (3). Two variables combined 
indicates a cross product. 




Table 3. Analysis of variance, regression coefficient’s, and statistics of fit for the 
dependent variable growth rate of McNair Vita-Graze Rye in Clemson greenhouse, 
January-.April 1974. 


Source 


df MS 

F 

Prob^F 

R* 

Regression 

Error 

Corrected total 

‘ - 

7 0.005707 
39 0.000769 
46 

7.42 

0.0001 

0.571 



Partial Regression 

Student's t 

Probability 


Coefficients 

for H0:B=O 

. >!H . 

Intercept 

0. 14169892 

6.53 

0.0001 

ST03/ 1000 

-0.01763732 

-2.88 

0.0065 

CUMRY 

0.00714499 

4.65 

0.0001 

DH3 

-0.00025372 

-2.07 

0.0451 

TN23/ 1000 

-0.00960439 

-2.82 

0.0075 

LN12 

-0.00000153 

-4.70 

0.0001 

TX3LN1 

0.00001761 

4.72 

0.0001 

TX3RH1 

0.00000148 

2.88 

0.0064 


The first number following each basic variable indicates the lag period. The second 
number indicates that the variable is squared (2) or cubed (3). Two variables combined 
indicates a cross product. 



Table 4. Analysts of variance, regressic* coefficients, and statistics of fit for the 
dependent variable growth rate of Blueboy wheat in Clemson field, November 1971- 
March 1972. 


Source 

df 

MS 

F 

Prob.> F 

R 2 

Regression 

8 

0.05875 

96.468 

0.0001 

0.859 

Error 

127 

0.00061 




Corrected total 

135 







Partial Regression 
Coefficients 

Student's t 
for HO: B=0 

Probability 
> 1*1 

Intercept 

0.05516383 

1.48 

0.1403 

TX12 

0.00020298 

11.05 

0.0001 

TN02 

0.00015460 

3.28 

0.0013 

TN22 

-0.00009994 

-1.83 

0.0696 

ESM22 

-0.00000442 

-1.22 

0.2254 

DH 23/1 000 

-0.00000292 

-3.45 

0.00C8 

DH 13/1000 

0.00000146 

1.82 

0.0710 

TN0RH3 

0.00000129 

2.95 

0.0037 

DH2RH3 

0.00000016 

2.82 

0.0056 


The first number following each basic variable indicates the lag period. The second 
number indicates that the variable is squared (2) or cubed (3). Two variables combined 
indicate a cross product. 



Table 5. Analysis of variance, regression coefficients, and statistics of fit for the 
dependent variable growth rate of Keowee barley in Clemson field, November 1971- 
March 1972. 


Source 

df 

MS 

F 

Prob.> F 

R 2 

Regression 

7 

0.05009 

44.9982 

0.0001 

0.854 

Error 

54 





Corrected total 

61 







Partial Regression 
Coefficients 

Student's t 
for HO: B = 0 

Probability 

> M 

Intercept 

-0.01871615 

-1.31 

0. 1956 

RH1TX1 

0.00000381 

4.35 

0.0001 

RH1TN0 

0.00000170 

1.79 

0.0784 

RH1DH1 

-0.00000113 

-3.20 

0.0023 

LN1ESM1 

-0.00007089 

-3.76 

0.0004 

LN1ESMQ 

0.00007193 

3.80 

0.0004 

TX1DH1 

0.00007634 

3.62 

0.0006 

TN0DH1 

0.00005901 

3.16 

0.0026 


The first number following each basic variable indicates the lag period. The second 
number indicates that the variable is squared (2) or cubed (3). Two variables combined 
indicate a cross product. 


Table 6, Analysis of variance, regression coefficients, and statistics of fit for the 
dependent variable growth rate of Wren's Abruzzi rye in Clemson field, November 
1971 -March 1972. 


Source 

df 

MS 

F 

Prob.> F 

R 2 

Regression 

8 

0.05946 

40.703 

0.0001 

0.723 

Error 

125 

0.00146 




Corrected total 

133 







Partial Regression 

Student's t 

Probability 


Coefficients 

for HO: B = 0 

>1*1 

Intercept 

-0.00164599 

-0.17 

0.8678 

TX1 

0.00314308 

3.26 

0.0015 

LN1 

0.00006702 

2.28 

0.0241 

TN32 

-0.00010073 

-2.10 

0.0377 

DH03/1000 

-0.00000802 

-6.40 

0.0001 

LN03/1000 

-0.00056888 

-3.12 

0.0023 

TN0TX1 

0.00032625 

5.30 

0.0001 

DH0TX2 

0.00001083 

1.68 

0.0945 

LN0TX2 

0.00000789 

3.02 

0.0030 


The first number following each basic variable indicates the tag period. The second 
number indicates that the variable is squared (2) or cubed (3). Two variables combined 
indicate a cross product. 



Table 7. Analysis of variance/ regression coefficients, and statistics of fit for the 
dependent variable, day of 50% planted in North Dakota, 1967-72, 1974. 


Source 

df 

MS 

F Prob> F 

R2 

Regression 

4 

1918.5994 

59.7372 0.0001 

0.8129 

Error 

55 

32.1173 



Corrected total 

59 






Partial regression 

Student's t 

Prob.>|t| 



coefficients 

for H0:B = 0 

Intercept 


136.6887 

54.99 

0.0001 

Product of E0 

N9S 

0.0055 

12.84 

0.0001 

Product of • v - C 

~nd X9S 

-0.1177 

-7.54 

0.0001 

Product of LOC and X3S 

0.2777 

6.24 

0.0001 

product of PP and X3S 

-0.0138 

-4.71 

0.0001 


N9S - Running nine day sum of average minimum temperature values (C°) 
X3S - Running three day sum of average maximum temperature values (C°) 
X9S - Running nine day sum of average maximum temperature values (C°) 
LOC - Location (north, central, south) 

PP - Preseason precipitation 
E0 - Estimated soil moisture 



Table 8.' Analysis of variance, regression coefficients, and statistics of fit for the 
dependent variable, day of emergence in North Dakota, 1967-72, 1974, 


Source 

df 

MS 

F 

Prob . > F 

R 2 

Regression 

8 

749. 19227 

37.495 

0.0001 

0.852 

Error 

52 

19.98114 




Corrected total 

60 







Partial regression 
coefficients 

Student's t 
for H0:B = 0 

Prob> |t| 

Intercept 

-1.026 X 10 2 

-1.435 

0. 1573 

E0 

9.257 X 10° 
-1.044 X 10" 1 

3.063 

0.0035 

E0 squared 

-2.516 

0.0150 

E0 cubed 

3.582 X 10" 4 
3.208 X 10' 1 

1.928 

0.0593 

N9S 

7.777 

0.0001 

PP squared 

8.364 X 10“ 2 

6.313 

0.0001 

Product of PP and N3S 

3.584 X 10‘ 2 

4.507 

0.0001 

Product of PP and X9S 

-1.564 X 10' 2 
-2.281 X 10" 1 

-7. 121 

0.0001 

Product of LOC and N3S 

-5.496 

0.0001 


EO - Estimated soil moisture 
PR - Preseason precipitation 
LOC - Location (north, central, south) 

N3S - Running three day sum of average minimum temperature values (C°) 
N9S - Running nine day sum of average minimum ter.iperature values (C°) 
X9S - Running nine day sum of average maximum temperature values (C°) 
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Fig. 15. Relationship of predicted date on which wheat emerged to actual date for 
the period 1967-72,. 1974, based on equation from Table 8. 







Abbreviations used in Appendix Tables 


GRO - Daily growth rate 
CUM - Cumulative growth rate 
FTXO “ Maximum temperature (°F ) 
FT NO - Minimum temperature (°F) 
LO - Langleys 

PO - Precipitation (inches) 

EO - Estimated soil moisture 


Appendix Table 1. Spring wheat growth rate and environmental factors at Williston, 
North Dakota, 1974 (ASCS field) 
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Appendix Table 3 . Spring wheat growth rate and environmental factors at Dicki 
North Dakota, 1974 
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Appendix Table 4. Spring wheat growth rate and environmental factors at Minot 
North Dakota , 1974 
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Appendix Table 5. Spring wheat growth rate and environmental factors at Clemson, 
South Carolina, 1974 
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